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Abstract 



We present experimental results for the ground state electrochemical po- 
tentials of a few electron semiconductor artificial molecule made by vertically 
coupling two quantum dots, in the intermediate coupling regime, in perpen- 
dicular and parallel magnetic fields up to B ~ 5 T. We perform a quantitative 
analysis based on local-spin density functional theory. The agreement between 
theoretical and experimental results is good, and the phase transitions are well 
reproduced. 
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I. INTRODUCTION 



Semiconductor quantum dots (QD's) are widely regarded as artificial atoms with proper- 
ties analogous to those of 'natural' atoms. Furthermore, systems composed of two QD's, 'ar- 
tificial' quantum molecules fQM's), coupled either laterally or vertically, have recently been 



vertically coupling two well defined and highly symmetric QD'sQ are ideally suited to study 
QM properties. We recently reported the addition energy spectra at zero magnetic field for 
such QM's function of interdot coupling strength.! 

In this work we present experimental and theoretical ground state electrochemical po- 
tentials for a diatomic QM in the intermediate coupling regime corresponding to an interdot 
distance, b = 3.2 nm, for magnetic fields (B) up to about 5 T. We assume here that the quan- 
tum mechanical coupling is sufficiently strong that the QM can be regarded as a symmetric 
'homonuclear' diatomic qmJ§ We consider two different configurations, one corresponding 
to an applied magnetic field parallel (B») to the drain current id flowing through the con- 
stituent QD's, and another corresponding to an applied magnetic field perpendicular (B±) 
to Id (see Fig. [I]). The latter has received relatively little attention.!'0Hl! We note that the 
QM physics we discuss in both magnetic field configurations is particularly relevant to the 
subject of solid-state quantum computing.0 

The interpretation of the experimental results here is based on the application of local- 
spin density-functional theory (LSDFT).i'E10 In the Bn case it follows the development of 
the method thoroughly described in Ref. |6], which includes finite thickness effects of the 
dots, and uses a relaxation method to solve the partial differential equations arising from a 
high order discretization of the Kohn-Sham (KS) equations on a spatial mesh in cylindrical 
coordinates (axial symmetry is assumed). To describe the less-common B± three- 
dimensional (3D) LSDFT code has been developed to handle configurations without any 
spatial symmetry. 

This work is organized as follows. In Sect. II we describe the experimental setup; in 
Sect. Ill we outline the method used to implement LSDFT in our axially symmetric QM 
system, and discuss the fully 3D configurations; and in Sec. IV we give the experimental and 
theoretical results. The interpretation of these results, and a short summary is presented in 
Sect. V. 



ically, two QD's which individually display clear atomic-like featuresJ3tl For the materials 
we typically use, the energy splitting between the bonding and anti-bonding sets of single 
particle (s.p.) molecular states, A S as, can be varied from about 4.5 meV (strong coupling) 
to about 0.1 meV (weak coupling) Bu In this paper, b, the thickness of the central barrier 
separating the two dots is fixed at 3.2 nm (A S as ~ 3 meV). Because this corresponds to 
intermediate coupling, we can reasonably neglect a small mismatch (of energy < Asas) be- 
tween the two dots, i.e. the QM is assumed to be symmetric ' homonuclear '.□ Figure [l] shows 
(a) a schematic section of a sub-micron circular mesa, diameter D, containing two vertically 
coupled QD's, and (b) a scanning electron micrograph of a typical mesa after gate metal 
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deposition. The starting material, a special triple barrier resonant tunneling structure, and 
the processing recipe are described elsewhere.01i^M Drain current 1^ flows through the two 
QD's between the substrate contact and grounded top contact in response to voltage Vd ap- 
plied to the substrate, and gate voltage V g on the single surrounding gate. The one structure 
we describe here (D ~ 0.5 microns) is cooled to about 300 mK or less. 



III. THEORY 

To analyze the experiments we have modeled the QM by two identical QD's stacked in 
the z direction (parallel to Id)- In this direction the QM is confined by two identical quantum 
wells of width w and depth V separated by distance b= 3.2 nm. We have taken V = 225 
meV and w — 12 nm, which are also appropriate for the actual experimental device. To 
improve on the convergence of the 3D calculations, the somewhat ideal sharp double-well 
profile has been slightly rounded-off, as shown in Fig. ^|. Given that the energy profiles of 
real structures are never abruptly sharp, this rounding-off is actually not unrealistic. 

In the remaining two directions perpendicular to z the QM is confined by a harmonic 
oscillator potential Vh = muj 2 (x 2 + y 2 )/2 of fixed strength huo = 4.42 meV. This lateral 
confinement energy has been determined for N = Q electrons using a general law0 that 
quantitatively describes the phases of QM's in the strong, intermediate and weak coupling 
regimes as a function of B\\ for a number of electrons, N, between 12 and 36. Lacking a 
better prescription at smaller ivJUl u has been kept fixed for all N analyzed here, and in this 
sense we present a parameter-free calculation. In the following we will denote by V c f(x, y, z) 
the total confining potential obtained by adding the double well profile to the harmonic 
oscillator potential, V^. We stress that V c f(x,y, z) is axially symmetric around the z axis. 

As it is well known, within LSDFT the ground state (g.s.) of the system is obtained by 
solving the Kohn-Sham (KS) equations. In the B\\ case the problem is greatly simplified by 
explicitly using the axial symmetry of the system. The additional terms in the KS equations 
due to the presence of an arbitrary magnetic field are given below. The inclusion of these 
terms crucially does not break the axial symmetry of the KS Hamiltonian in the Bn case. 

In the symmetric gauge the vector potential A(f) corresponding to a constant magnetic 
field B is written as A = (B A f )/2, and its contribution to the KS Hamiltonian is 

H m = 7 ^B-L + -^A 2 + g ; l i B B-S , (1) 
2mc 2mc z 

where g* is the effective gyromagnetic factor, L and S respectively are the orbital and spin 
angular momentum operators, and /ie is the Bohr magneton. Writing B = B(sm0B, 0, cos#b) 
and introducing the cyclotron frequency u c = eB/mc, it can be easily checked that Ti m = 
T-L mR +iHmj, with 
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where rj a =+l(— 1) for er=T(!) with respect of the direction of the applied magnetic field. 



We have used effective atomic units % = e 2 /e 



m 



-1, where e is the dielectric constant, 



and m the electron effective mass. In units of the bare electron mass m e , m 



m m P 



In 



a e/m* with a = % /m e 



this system, the length unit is the effective Bohr radius Oq 
and the energy unit is the effective Hartree H* = Hm*/e 2 
the following values: g* = —0.44, e = 12.4, and m* = 0.067. This yields 
H* ~ 11.86 meV. From now on we will write the equations in these units. 

Equation (0) reduces to the B» case when Q B — 0, and to the B± case when 6 B 



For a QD in GaAs, we take 
97.94 A and 

7r/2. In 



the former, since (x-j^ ~U~§^) xs proportional to L z , the problem remains axially symmetric. 
A detailed description of how the KS equations have been solved in this geometry can be 
found in Ref. ||. 

In 3D the KS equations read 



d 2 



D 2 



d 2 



2 \dx 2 dy 2 dz 2 



+ Vcf{x,y,z) 



(3) 



+ V H + V xc + W xc Va + H m ] y, z) = e a V a (x, y, z 



The expression in the brackets is the KS Hamiltonian 'Hks-, an d V H (x, y,z) is the direct 
Coulomb potential. V xc = 5£ xc (n, m)/5n\ g s ^ and W xc = 5£ xc (n, m)/5m\ g s , are, respectively, 
the variation of the exchange-correlation energy density S xc (n, m) in terms of the electron g.s. 
density n(x,y,z), and of the local spin magnetization m(x,y,z) = n^x^jz) — n^x^.z). 
The exchange-correlation energy has been taken from Perdew and Zunger,li3 and £ xc (n,m) 
has been constructed as indicated in Ref. |6|. It is worth noticing that if B ^ then the s.p. 
wave functions ^ a (x, y, z) are complex, with their real and imaginary parts being coupled 
by H m . 

The KS and Poisson equations are solved on a 3D mesh after discretizing them using 
7-point formulas, and using a two-grid version of the one-way multigrid method described in 
Ref. [19]. The Poisson equation is solved using a first order relaxation schemed The required 



value of the Coulomb potential at the mesh boundary is obtained by a standard multipole 
expansion up to eighth order. The KS equations are solved using an imaginary time method, 
involvingthe third-order expansion of the forward solution of the imaginary time diffusion 
equationcj 



-{n KS - e)tt 



(4) 



i.e., 



*(r + 8t) - *(t) = A*(r) = -5t(H K s - e)*(r) 



2 6 



(5) 



where e = (^(t)\Hks\^( t )) ■ To further accelerate the self-consistent solution of both the 
KS and Poisson equations, we use the preconditioning smoothing operation described in Ref. 
21. In the KS case, this means that A\l/(r) has been smoothed as proposed in this reference. 
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The perfomance of the code has been further improved by adding a 'viscosity term', i.e., Eq. 
(|) has been changed into 

*(r + 8t) - *(r) = A*(r) + a v [*(r) - *(r - <5r)] . (6) 

A viscosity term has also been included in the solution of the Poisson equation. We have 
used a 45 x 45 x 67 mesh with spatial steps Ax = Ay = 5.67 nm, and Az = 0.89 nm. The 
large asymmetry between the spatial meshes is motivated by the sharpness of the confining 
potential in the z direction. The heuristic viscosity parameter ay is fixed to a value of 0.8, 
and the time step 5t to the value of (Az) 2 in effective atomic units. The stability of our 
results against the increase of the number of mesh points has been checked, and at B = 
we have used the results of the axially symmetric code to test the results obtained with the 
3D code. 

The accuracy of LSDFT for the B values of interest has been assessed by comparing the 
results for a single QD with those obtained using the current spin-density functional theory 
(CSDFT), which in principle is better suited to high magnetic fields than LSDFT. Since 
CSDFT is a two dimensional (2D) theory,!! we have also compared our LSDFT results with 
those obtained using the 2D-LSDFT which is implicit in any implementation of CSDFT, 
in particular see that of Ref. |2^. The low and high field borders of the maximum- density 
droplet (MDD) phase using strictly 2D-LSDFT and CSDFT, as described in Ref. ^3], have 
been obtained for QD's with N = 20, 28 and 36 electrons laterally confined by a harmonic 
oscillator potential of energy fvuj = 7.6 iV -1 / 4 (meV). Such parameterization of the confining 
potential within LSDFT reproduces the experimental MDD.0 The results are shown in 
Table |, together with the values obtained by using the present 3D-LSDFT for the same 
lateral confining potential. From Table | we can see that the overall agreement between the 
three calculations is clearly good, and thus we can confidently use LSDFT in the present 
calculations. 



IV. RESULTS AND DISCUSSION 



The experimental ground state electrochemical potentials for N = 3 to 6, as a function 
of B, are shown in Fig. |3| for (a) the parallel case, and (b) the perpendicular case. What 
is actually shown is the S-field dependence of the 3rd, 4th, 5th and 6th current (Coulomb 
oscillation) peaks measured by sweeping V g in the linear conductance regime for a small 
Vd — 0.1 mV. It is clearly evident that the dependencies for parallel and perpendicular cases 
are very different - in particular the former is stronger than the latter. We now attempt to 
explain the general appearance in both cases, and in particular the features marked by the 
different symbols, by using the computational methods described in the previous section. 

We first address the B\\ case. As in single QD's, at low £>-fields, upward kinks (cusps) 
in the experimental N electron g.s. QM electrochemical potentials as a function of B are 
interpreted as changes in the N electron g.s. configuration of the QM which arise from s.p. 
level crossings. i@M We have plotted in Fig. |4] the calculated g.s. electrochemical potential 
fJ>{N), defined as 

/i(A) = U(N) - U(N - 1) , (7) 
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where U is the total energy of the N electron QM g.s., as a function of Bu up to A^ = 6. 
To label the g.s. configurations we have used the usual notation of molecular physics for 
s.p. electronic orbitals.a Upper case Greek letters are used for the total orbital angular 
momentum. We have also used the adapted versio ni of the ordinary spectroscopic notation 
+1 Lg U with S being the total \S Z \, and L being the total \L Z \. The superscript +(— ) refers 
to even (odd) states under reflection with respect to the z = plane bisecting the QM. Even 
states are bonding (symmetric) states, and odd states are anti-bonding (anti-symmetric) 
states. The subscript g(u) refers to positive (negative) parity states. All these are good 
quantum numbers in the B\\ case and can be used to label the different g.s.'s ('phases'). 
Following Refs. we have also calculated the 'isospin' quantum number (the 'bonding 
number' in molecular physics) defined as I z = (Nb — Nab)/2, with Nb(ab) being the number 
of electrons in bonding (anti-bonding) s.p. states. This is an exact quantum number for 
homonuclear QM's in a parallel magnetic field.0 

Given the complexity of real vertical QM structures and the challenge in modelling 
them, a comparison between Figs. |3|(a) and |] reveals a rather good agreement between 
theory and experiment. As a guide, and consistent with the calculated states and the 
observed Bu dependence, we indicate in Fig. |3](a) in simple box styled the dominant g.s. 
configurations at or near B = 0, and others at higher field which are stable over a relatively 
wide range of Bn . Up and down arrows indicate spin- up and spin-down electrons, and black 
(grey) arrows represent electrons in bonding (anti-bonding) s.p. states@@ For N = 3, 5, 
6, near B = 0, because the g.s.'s are close to each other, i.e. stable over a fairly narrow 
range, we show two configurations which in practice are hard to resolve. Some of these 
involve the population of the lowest anti-bonding state with a single electron, so isospin 
is non-maximal. Above B = 1 T, however, all the anti-bonding states are depopulated so 
isospin is maximal (I z = N/2), and filling of the QM resembles that of a single QD. The 
identifiable g.s. transitions in Figs. §(a) are marked by black triangles. As expected, most 
appear as upward kinks. A couple, see first kinks for N = 5 and 6, appear as "downward" 
kinks because of the g.s. transitions which occur at almost the same Bu in A^ = 4 and 5 
respectively. 

Looking further at other details in Fig. |], for N = 2, the singlet-triplet transition 
occurs at about 4.6 T which is close to the experimental valueEl of ~ 4.2 T (not shown). 
We have found from the calculations an MDD configuration made of electrons filling just 
bonding s.p. states (MDDb), which has a total angular momentum L z = N(N — l)/2, 
and extends from ~ 4.9 to ~ 9.5 T for N = 3, from ~ 5.1 to ~ 9.0 T for N = 4, from 
~ 5.4 to ~ 8.8 T for A^ = 5, and from ~ 5.6 to ~ 8.3 T for N = 6. These results are 
at variance with those of Ref. where an MDDb g-s. was found for N = 3, but not for 
larger values of N. The reason of this discrepancy may be attributed either to the strictly 
2D model used in their calculation to represent the constituent QD's, or more likely to 
their particular implementation of CSDFT.0 A quantitative description of the spin phases 
of QM's in different coupling regimes, along with the experimental data, will be presented 
elsewhere .0 

In Fig. |] we can see that for the larger N values studied here, the increase in angular 
momentum of the QM g.s. as it evolves from B = towards the MDDb is accompanied 
by two isospin flipsE^ caused by electrons jumping from anti-bonding to bonding states and 
vice-versa. Phase transitions from — to + g.s.'s involve AI Z = 1 flips, whereas those from 
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+ to — g.s.'s involve AI Z = — 1 flips, and both are clearly seen in Fig. |4] for N = 5 and 6. 
Interestingly, they only happen for B\\ < 2 T. We can see that after reaching the i/g = 2 g.s. 
(i.e., a filling factor two QM state made of just bonding s.p. states), which corresponds to 
the X A+ phase for N = 4, to the 2 r+ phase for N = 5, and to the 1 J+ phase for N = 6, only 
bonding s.p. states are occupied, and as a consequence the QM reaches the MDDb state in 
a similar way to how a single QD reaches the MDD state, namely by populating bonding 
s.p. states of higher and higher s.p. orbital angular momentum I values.^ In general, these 
isospin flips can produce a complex pattern in the s.p. spectrum as a function of B». As an 
example of this complexity, we present in Fig. [5] the s.p. levels for N = 6 as a function of I 
for different B» values. It can be seen in this figure that as B increases, the QM undergoes 
isospin flips. Firstly, the / = 0| anti-bonding s.p. state becomes occupied, as shown in 
the panels corresponding to B = 0.5 and 1.2 T. After another isospin flip caused by the 
depopulation of the same s.p. state, the QM reaches the v B = 2 phase corresponding to 
the configuration {B = 3 T panel). From this phase on, the spin polarization steadily 
increases until the QM reaches the MDD# phase (B = 6 T panel). 

In the B± case, crucially, within LSDFT, the only good quantum number is the spin 
projection along the direction of the applied magnetic field, which we call s±, and the g.s. 
electrochemical potentials as a function of B± are expected to be much smoother than in 
the Bu case. This situationHEiHl! unlike the Bn case, lacks an analytical solution even for 
the case of non-interacting electrons. We show in Fig. |]the calculated non-interacting s.p. 
spectrum as a function of B±. At B = the energy difference between bonding and anti- 
bonding I = s.p. states is just As as (likewise for the / = 1 states). The energy difference 
between I = 1 and bonding (or anti-bonding) states is just Hcu. Similar results have been 
reported elsewhere. Note that the effective value of A S as is reduced as B± is increased 
because the tunneling electron trajectories are increased in the barriers. The small splitting 
between spin up and down states that originate from a common s.p. state with well defined 
orbital angular momentum at B = are due to the Zeeman term. Increasing B± has also 
an effect on the quantum mechanical coupling of the QD's forming the QM, as the already 
sizeable electron density in the middle of the central barrier increases further (see n(z) in 
Fig. |2|) which magnifies exchange-correlation effects between both QD's. The compression 
of n(z) in each well as B± increases shown in Fig. || can be understood from Eq. (H): note 
that when 6b = 7r/2, the first term in 7im R becomes u 2 (y 2 + z 2 )/8. This also produces an 
upward shift of the s.p. spectrum of the QM as a whole (the horizontal dashed line in Fig. 
|2] lies above the horizontal solid line in the same figure). 

The even (bonding) or odd (anti-bonding) character of the s.p. levels defining a QM state 
is strictly lost when a magnetic field perpendicular to is present. Intriguingly, however, 
the bonding/ anti-bonding character present at B = is sometimes retained to a large degree 
by the s.p. states at finite B± values. We have indeed found that the expectation value of 
the z —>■ — z reflection operator 

(IT) = J df^f*(f)IL z ^f(f) = J df^/*(x,y,z)^(x,y,-z) (8) 

is very close to ±1, as it should be for bona fide bonding/anti-bonding states, in many cases 
even for relatively large values of Bj_. 

As an example of this, we show in Fig. |7| the energies of the occupied s.p. states as 
a function of for N = 5 and 6. Solid triangles represent 'quasi-bonding' states with 
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{U z ) > 0.95. Note that at T for iV = 5 the s.p. bonding state at e ~ 48.2 meV is two 
fold degenerate, and likewise for two of the N = 6 s.p. bonding states at e ~ 52 meV. Open 
triangles represent 'quasi-anti-bonding' states with (IT 2 ) < —0.95. Actually, there is only 
one such occupied anti-bonding s.p. state for N = 5 at B = 0, and none for N = 6. All 
other open symbols (circles and squares) correspond to s.p. states with negative (II Z ) values 
larger than -0.95, i.e., cannot really be regarded even as 'quasi-anti-bonding' states. 

The figure also shows that states that evolve from I = s.p. states at B = retain a 
quasi-bonding character up to quite high values of B± (at least up to 5 T), whereas other 
states, that at B = are I = 1 s.p. states, do not. The quasi-bonding robustness of the 
lower-lying s.p. states may be due to the small effect that the applied magnetic field has 
on states that are I = s.p. states at B = 0. The B± evolution of what at T are the 
2p-states is rather similar for N = 5 and 6 with a change from solid to open symbols near 
Bj_ = 4 T. Interestingly, in spite of the lack of any spatial symmetry in the system when 
a perpendicular field is applied, the s.p. levels are still clearly distributed into shells as in 
the non-interacting caseBO Notice also the different splitting between f and j states. For 
saturated (zero) spin (N = 6 case), this is essentially due to the small Zeeman term, whereas 
for non-saturated spin (N = 5 case) the splitting is mostly dueil to the spin- dependent part 
of the exchange-correlation energy, W xc term in Eq. (^), and this effect is larger the higher 
the value of the g.s. spin. This explains the sizeable splitting between the two lower lying 
s.p. levels for N = 5 up to B± ~ 0.5 T and the splitting of all the s.p. levels for N = 6 
above B± ~ 3.5 T (see also Fig. ||). 

The calculated g.s. electrochemical potentials are shown in Fig. § as a function of Bj_. 
Comparing with Fig. 0(b) it can be seen that the agreement with experiment is good for 
3 < N < 6. We have indicated the value of the total S± for all the relevant g.s. phases. In 
the B = to 5 T range, there are some B± induced changes in S±, and these give rise to 
upward kinks (also marked in Fig. 0(b) by solid down triangles). Some downward kinks, 
identified by vertical arrows in the N = 5 and 6 g.s. electrochemical potentials do not 
correspond to changes in the A-electron S±. They are associated with s.p. level crossings 
between s.p. states of different (U z ) value of the N — 1 electron system. This is the case for 
N = 5 at B± ~ 1.5 T and N = 6 at Bj_ ~ 0.75 T, as can be seen in Fig. [8] (also marked 
in Fig. |3](b) by solid circles). Because of the lack of spatial symmetry in the system, we do 
not, in general, attempt to identify the (dominant) g.s. configurations. The configurations 
shown in simple box style in Fig. 3 (b) are the dominant g.s. configurations at T and they 
are expected to remain so for small values of Bj_. The singlet-triplet transition for N = 2 
in Fig. 8 appears at ~ 4.7 T, a value compatible with that found in the B» case. The Bj_ 
induced sing let-triplet transition in the experimental data is discussed elsewhere.!! It can 
also be seen that for N = A Hund's first rule like filling occurs for B± < 1.5 T, even if the 
g.s. configuration is not strictly axially symmetric. Nonetheless, for iV = 6 at B = we 
have found an axially symmetric configuration corresponding to a 2D harmonic oscillator 
shell like filling. 

V. SUMMARY 

We have thoroughly discussed the ground state electrochemical potentials of a few- 
electron semiconductor artificial QM in the intermediate coupling regime. A detailed com- 
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parison between experimental data and LSDFT calculations shows overall a good agreement 
for both parallel and perpendicular magnetic fields. The agreement is even more remarkable 
since the frequency u of the lateral confining potential has not been used here as a fitting 
parameter, but rather it has been derived from a law strictly valid for larger values of N. 
Had we used an even smaller value of u, the agreement would have been even better. 

Any sensible comparison with the results of other calculations and with the experi- 
mental data should consider the strong influence of the confinement on the actual g.s. QM 
phases. Large u values obviously favor the occupation of anti-bonding states. Consequently, 
decreasing uj might 'wash-out' phase transitions involving anti-bonding states. For example, 
see the 2 n+ — > 2 A~ transition for N = 5 or the x £+ — > 3 U g transition for N = 6 in Fig 
4. We have checked that this is indeed the case when the N dependent confining potential 
htu = 5.78 A _1//4 (meV) of Refs. |3|,|26| is used. The same is true for the 2 X~ — > 2 n+ transition 
for iV = 3 (these two states are practically degenerate at B = 0). Thus, a comparison be- 
tween theory and experiment may help to find accurate and realistic values for the effective 
lateral confining potential. 

Analysis of the B± case has shown that the N = 2 singlet-triplet transition sensibly 
occurs at a similar B value to that in the Bn case. In spite of the absence of the strong 
Landau quantization inherent to the Bn situation, the s.p. levels are, to a large extent, still 
distributed into shells. We have also found that spin effects arising from the spin dependence 
of the exchange-correlation energy can dominate those caused by the small Zeeman term as 
is also the case for a single QD. 
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29 The simple schematic box notation in Fig. |3] that represents the molecular states is more 
intuitive than the full spectroscopic notation used in Fig. ||, although we do not use it 
in Fig. |] because of lack of space. Nonetheless, it is easy to go from one notation to the 
other since each box represents a s.p. level characterized by a radial quantum number n 
and an angular momentum I. The black (grey) arrows denote the thrid component of spin 
for electrons in bonding (antibonding) states. At T, for both bonding and antibonding 
states separately, the lowest box is the (n, I) = (0, 0) Is state, and above that are two 
boxes, degenerate in energy, which are the (0, ±1) 2p states (see for example N = 5 in 
Fig. |3](a)). For the spectroscopic notation, the superscript — of states shown in Fig. |] is 
associated with an electron in an antibonding box in Fig. 0(a), and the total S z and L z are 
readily obtained by adding the spin and orbital angular momentum of all the electrons. 
Parity (g or u) is obtained by multiplying the parities of the s.p. states. 

30 We want to stress that in the B\\ case we have explicitly imposed the s.p. wave functions 
to be axially symmetric, but have not forced them to be specularly symmetry about the 
z = plane. Interestingly, in the case of a two-electron homonuclear QM in the weak 
coupling limit corresponding to fairly large interdot distances, we have foundi that the 
LSDFT g.s. is a broken-symmetry state made of two s.p. wave functions each of them 
localized in one of the constituent dots, in a singlet spin state. The specularly symmetric 
configuration has a slightly higher energy and corresponds to a triplet spin state. 

31 The authors of Ref. |26| have identified the vorticity V with the external magnetic field 
B. While this is indeed the case for a homogeneous 2D electron gas, and a fairly good 
approximation for large QD's, this may not be a good approximation for small systems, as 
in the present case. CSDFT implementations better suited for small systems are presented 
in Refs. 22-25, for example. 

32 Due to the convention adopted for the s.p. wave functions,! the s.p. orbital angular mo- 
mentum is actually —I. 
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TABLES 



TABLE I. Comparison between two-dimensional CSDFT and LSDFT, and three-dimensional 
LSDFT results for one single QD. L(R) denotes the left(right) border of the MDD phase in the 
N — B\\ plane. 





2D-CSDFT 


2D-LSDFT 


3D-LSDFT 


N 


L(T) 


R(T) 


L(T) 


R(T) 


L (T) R (T) 


20 


5.6 


6.3 


5.4 


6.5 


5.4 6.4 


28 


5.6 


6.1 


5.5 


6.2 


5.5 6.1 


36 


5.6 


5.9 


5.5 


6.0 


5.6 6.0 
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FIGURES 



FIG. 1. (a) Schematic diagram of mesa containing two vertically coupled quantum dots and 
(b) scanning electron micrograph of a typical circular mesa. 

FIG. 2. Double quantum well potential used in the calculations. Two electronic densities n(z) 
corresponding to the N = 6 QM for B = (solid line) and Bj_ = 5 T (dashed line) are shown. 
n(z) has been obtained by integrating the electronic density n(x, y, z) over the x and y coordinates. 
The energy of the corresponding occupied upper lying s.p. level is also represented by a horizontal 
solid line for B = 0, and by a horizontal dashed line for B± = 5 T. 

FIG. 3. Experimental f?-field dependence of the 3rd to 6th Coulomb oscillation peaks (g.s. 
electrochemical potentials for 3 < N < 6) in (a) B\\ case, and (b) Bj_ case. 

FIG. 4. Theoretical ground state electrochemical potentials in the Bu case (dots) for N < 6. 
The lines have been drawn as a to guide the eye. The vertical ticks along the fi(N) lines indicate 
phase boundaries. The various states are identified by standard spectroscopic notation discussed 
in the text. 

FIG. 5. Single particle energy levels as a function of I for different values of Bn at 
N = 6. Upward (downward) triangles denote |(|) spin states. Open(solid) triangles correspond to 
anti-bonding(bonding) states. The horizontal lines represent the Fermi level. The value of -By is 
indicated in each panel. 

FIG. 6. Energies of the nine lower-lying non-interacting s.p. levels as a function of B±. As as 
and too are marked. For each symbol, the direction of is indicated in the box. 

FIG. 7. Single particle energy levels as a function of B± for N = 5 (left panel) and 6 (right 
panel). For each symbol, the direction of s± is indicated in the box. States two fold degenerate 
are indicated by x2 symbol. Solid and open symbols are discussed in the text. 

FIG. 8. Theoretical ground state electrochemical potentials in the B± case (dots) for N < 6. 
The lines have been drawn as a to guide the eye. The vertical marks along the fJ.(N) lines indicate 
phase boundaries. The value of S± in each phase is given. We have indicated by vertical arrows 
downward kinks arising from s.p. level crossings in the iV-1 electron ground state that do not 
produce phase transitions in the N electron ground state. 
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